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Abstract 

The new light-front coupled-cluster (LFCC) method for the nonperturbative solution of Hamil- 
tonian eigenvalue problems is described and then illustrated in an application to quantum electro- 
dynamics. The method eliminates any necessity for a Fock-space truncation and thereby avoids 
complications associated with such a truncation. An LFCC calculation of the electron's anomalous 
magnetic moment is formulated for a truncation that, for simplicity, excludes positrons but keeps 
arbitrarily many photons. 
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I. INTRODUCTION 

In order to have a better understanding of hadronic physics within quantum chromo- 
dynamics, it is imperative to have methods for the nonperturbative solution of the theory. 
Lattice gauge theory [1] has been quite successful in this regard. The use of Dyson-Schwinger 
equations [2| has also born fruit. Both, however, are somewhat indirect, in that they do not 
provide Minkowski-space wave functions from which one could compute observables directly. 
An alternative that does deal directly with such wave functions is the light-front Hamilto- 
nian method [3j. There wave functions appear as coefficients in Fock-state expansions of the 
eigenstates. 

The light-front Hamiltonian approach now has long history [3|, beginning with Dirac's 
initial formulation of the coordinate choice [4|. A key ingredient in practical calculations 
has been the truncation of Fock space. This yields a finite set of equations for a finite set of 
wave functions, which are usually solved numerically. 

Unfortunately, such truncations cause several difficulties, including uncanceled diver- 
gences [5( , broken Ward identities [6| , and loss of Lorentz covariance [7|, |8j and gauge imparl- 
ance [9] . These come from the nonperturbative analog of decomposing a Feynman diagram 
into time-ordered contributions and throwing away those contributions that involve more 
than some fixed number of intermediate particles. Some of these difficulties can be some- 
what mitigated by use of a sector-dependent parameterization 1CH12J]. but they are not 



eliminated and can, in fact, lead to ill-defined wave functions [13 . 

To avoid these difficulties, we have recently proposed a light-front coupled-cluster (LFCC) 
method [l4J . Fock space is not truncated. Instead, an eigenstate is constructed from a va- 
lence state |0) by the action of an exponentiated operator T. With inclusion of a normalizing 
factor VZ, the eigenstate is written as yZe T \(f)). The problem is then to find |0) and T. 
To make the problem finite, T is truncated, but the exponential of T is not, so that e T \4>) 
includes contributions from all relevant Fock sectors. The method relies on light-front coordi- 
nates, to make Fock-state expansions well-defined and to facilitate separation of internal and 
external momenta, and on the mathematical techniques of the many-body coupled-cluster 
method 15|, \l& , hence the name. 



The original coupled-cluster method is used primarily to solve many-body problems where 
the number of particles is large and unchanging. It was first developed by Coester and 



Kummel [15( for applications to the many-body Schrodinger equation in nuclear physics. It 
was extended to many-electron problems in molecules by Cizek [17|, which was eventually 
followed by extensive development in quantum chemistry [16|. The basic idea is to form an 
eigenstate as e T |0), where \<p) is a product of single-particle states, as in the Hartree-Fock 
approximation, and the terms in T annihilate states in |0) and create excited states, to build 
in correlations. A finite numerical calculation is then done by truncating T at some (small) 
number of excitations. There has also been a series of applications to field theory [18]. These 
rely on Fock-state expansions in equal time and necessarily focus on the structure of the 
vacuum. Particle states are built on this vacuum. There was some success in analyzing </> 4 
theory in 1 + 1 dimensions. 

Our purpose here is to illustrate the LFCC method in a concrete calculation in a gauge 
theory. We consider the state of the dressed electron and its anomalous magnetic moment in 
an arbitrary covariant gauge. The theory is regulated by the inclusion of Pauli-Villars (PV) 



photons and electrons [5|, |19J. We do truncate Fock space but only to exclude positrons; all 



possible single-electron, multi-photon states are retained. The exclusion of positrons is not 



required by the method, but is done to simplify the illustration. 

We begin in Sec. |TT]with an extended description of the LFCC method 14|. The appli- 
cation to QED is developed in the following two sections; Sec. IIHI focuses on the eigenvalue 
problem, and Sec. II VI on the calculation of the anomalous moment. The details of light-front 
QED in an arbitrary covariant gauge [9[ are discussed in Appendix |AJ including a new anal- 
ysis of the gauge projection. Details of the construction of the effective LFCC Hamiltonian 
are given in Appendix [B] 



II. A LIGHT-FRONT COUPLED-CLUSTER METHOD 

We use light- front coordinates [3, |4|, where time evolves along x + = t + z. The spatial 
coordinates are x = (x~,x±), with x~ = t — z and x± = (x,y). The light-front energy is 
p~ = E — p z and momentum p = (p + ,p±), with p + = E + p z and p± = (p x ,Py)- The mass- 
shell condition p 2 = m? becomes p~ = (m? + p 2 L )/p + - In a system with total momentum 
P, a constituent with momentum p is defined to have a longitudinal momentum fraction 

y = p + / P + and relative transverse momentum k± = p± — yP±. A creation operator <v(jp) 

for such a particle will be written as a^(y,kj_;P_). Clearly, we have a^(l,0;P) = a^(P). 
Given a light-front Hamiltonian V~ , we wish to solve the fundamental eigenvalue problem 

/If 2 , p2 

v-mp))= p + - mp)) (2.i) 

for the eigenmass M and the eigenstate \ip{P_)), in a basis where the momentum V_ is diagonal 
with eigenvalue P. The normalization is chosen to be 

(V>(p'M£)> = S(P! - P). (2.2) 

The eigenstate can be generically expanded in a Fock basis {\n;pi)} as 



mp)) = J2 II d wdSu.)*(i-X)w)*(Z)^-)^(w»Sx)l n ;a>' ( 2 - 3 ) 



where n is the number of constituents and p. their light-front momenta. The eigenvalue 
problem then becomes a coupled set of integral equations for the wave functions ip n . In 
this form, the problem requires truncation to a finite set of wave functions, and the various 
complications follow, as discussed in the Introduction. 

To avoid truncation of the Fock space, we construct the eigenstate as 

|V(P)} = yfZe T \cj>(P)). (2.4) 

The valence state |0(P)) contains a small number of constituents and will be obtained from 
diagonalization of an effective Hamiltonian in the corresponding valence sector, where the 
small number of particles is fixed. The operator T is constructed to preserve all the quantum 
numbers of the valence state, including momentum P and angular momentum projection 
J z , and to always increase the number of constituents. The exponential of T then generates 
all the higher Fock states and their associated wave functions. The normalizing factor yZ 
is chosen to allow the valence state to have the same momentum normalization as the full 

(<K£'M£)} = <*(£'-£)• (2-5) 



The eigenvalue problem can now be written as 

/l/f 2 _|_ p2 

V-\HP)) = p + ± \m)), (2-6) 

with V~ = e~ T- P~e T the effective Hamiltonian. By introducing a projection P v onto the 
valence sector, defined as the sector of Fock space containing the Fock states found in the 
valence state, we separate the fundamental eigenvalue problem into a valence eigenvalue 
problem 

/If 2 _i_ pi 

P V V-\<P(P))= p + + ^ |0(£)) (2.7) 

to be solved for the valence state, and a set of auxiliary equations 

(l-P v )V~\<P(P)) = (2.8) 

that determine the operator T. The equations must be solved simultaneously, of course. 

For an exact solution, this system of equations will not be finite; the T operator will have 
an infinite number of terms, and there will be a correspondingly infinite number of auxiliary 
equations to fix them. We therefore truncate T to a few terms and truncate the projection 
1 — P v to include only enough higher Fock states to have enough equations to solve for the 
terms in T. The construction of these equations is aided by use of the Baker-Hausdorff 
expansion 

W = V- + [P-,T] + ^[[V-,T],T] + ■ ■ • , (2.9) 

which can be truncated to the finite number of terms relevant for the truncated projection. 
The restriction that T always increase particle number is critical for this truncation of 
the Baker-Hausdorff expansion. The exponential of T is not truncated, so that e T |0) still 
contains all the higher Fock states. 

The infinity of Fock states included in e T |0(P)) makes the normalization step nontrivial. 
Direct computation of Z requires an infinite sum. However, we can still compute expectation 



values of observables, using a trick borrowed from the original coupled-cluster method [16 . 
Let O be the Hermitian operator representing the observable of interest. We wish then to 
compute 

(0) = Z(<p(P)\e Tt de T \ ( p(P)). (2.10) 

To do this, we define 

= e- T 6e T (2.11) 

and 

mE)) = Ze T 'e T \4>(P)) = VZe T '\^P), (2.12) 

so that 

(O) = (V>(£)|O|0(P)}. (2.13) 

By construction, we have 

<<K£'M£)> = mP!M(P)} = KEl-P) (2.14) 

and 

V^$(P)) = e T 'v-e- Tt ^Ze Tt \ij(P)) = 4Zt TX V'\^P)) = M p ^ P± [fe))- (2.15) 



Thus, we can find \ip(E)) by solving the eigenvalue problem for V ^ V , at the same mass 
M, and normalizing \i/)(P)) to the valence state \<f){E)) according to Eq. (I2.14p . When T is 



truncated, \ip(E)) must also be truncated [14]], to the valence sector plus those Fock states 



associated with T\(f>(P)). This leaves a finite set of coupled equations for the functions in 
\4>{E))- Also, the Baker-Hausdorff expansion of the effective operator O can be truncated 
at a finite number of terms. 

One useful generalization of this technique is the calculation of off-diagonal matrix ele- 
ments 

(MP2)\0\MPi)) = V%Z~2(MP.2)\e T *de Tl \MPi)- 

This can be done by defining Oi = e~ Ti Oe Ti and \ipi(E)) = Z^i e Ti |0j(P)) and then con- 
sidering the two rearrangements for an Hermitian O 

Z l ~ /^ \ — -T 2 7i| 



and 




(MP2)\0\MEi)) = \l^(MP2)\02e-^e^\MPi)) 

^2 



(UPi)\6\UP.2)) = MiMEi^e-^lM^)}- 

V ^1 

The unknown normalization factors cancel in the product of the first by the conjugate of 
the second, which yields 



(ME2)\0\MPi)} = ^J (MP2)\02e- T ^\MEi)} (2-16) 



xV(^i(£i)|Oie-^|0 2 GP 2 )>*. 

The correct phase can be obtained by checking one or the other of the individual matrix 
elements. The factors of e~ Tl e T2 and e~ T2 e Tl can be evaluated with use of power series 
expansions for the exponentials or the Zassenhaus expansion |20j 

gTa-Ti _ e 2b e -Ti e |pZi,2\] e i[31 1 [2],,2\]]+|[[21,,ri],Ib] _ _ _ ( 2 .17) 

Only a finite number of terms will contribute, because the 7] only increase particle number 
and the initial and final states include only a finite set of Fock sectors. In the special case 
of an expectation value, the formula in (12.161) reduces to (I2.13p . 

Another generalization of the technique for expectation values is to include a projection P s 
of the eigenstate onto a subspace. We have in mind applications to gauge theories where the 
projection is onto a physical subspace of states that satisfy the gauge condition; however, in 
this Section we will not need to be specific. Let \ip s (E)) — \fZ~sP s z T \ < t>(.E)) be the projected 
state, normalized such that 

(UE!)\UE))=5{E-E). (2.18) 

The unprojected state is related by 

l^(£)} = MpME))- (2-19) 



The expectation value of an operator is given by 

(0) s = (ME)\d\MP)) = ZME)\e Tt P}OP s e T \^P)). (2.20) 

On introduction of the unprojected left-hand eigenstate \ip(P_)), this becomes 

(0) s = §M£)|e- r P|OP s e T |0(P)>. (2.21) 

To obtain the ratio Z s /Z, we assume Z s and Z to be independent of the total momentum 
and use 

5{E-P) = (MP!M(P)) = ^(0(£')k T1 [ P}P s e T \m)) = ^${P!)\e- T PtP s e T \<l>{P))- 

(2.22) 
Integration over P' then yields 



|dP'(^(P')|e- r P]P s e T |0(P)). (2.23) 



Z 

We therefore obtain 

Q (^(P)|e^P;OP s e r |0(P)) 

/ ' dP!{iP{E)\e- T PlP s e T \<p(P))' 

The remaining matrix elements are computed from the known \<p(P_)) and \tp(P)) by first 
expanding the exponentials, keeping the total number of T factors consistent with the trun- 
cation of T, carrying out any contractions of annihilation and creation operators, and then 
applying the projections P s and Pj\ A Baker-Hausdorff expansion is not useful here, be- 
cause of the projections, but the truncated expansion of the exponentials is equivalent to a 
truncated Baker-Hausdorff expansion. 

III. AN APPLICATION TO QED 

To show how the LFCC method works for a gauge theory, we consider the dressed- 
electron state in QED, truncated to exclude positrons. The theory is regulated by one Pauli- 
Villars (PV) electron, with coupling coefficient (3% — 1, and two PV photons, with coupling 
coefficients £i and £2- Additional PV fields are not needed if positrons are absent J2l|. The 
construction of the covariant-gauge light-front Hamiltonian is discussed in Appendix |XJ with 
additional details given in J9|. The valence state is 

i^ ± (£)} = E^ ± A t ± (aio), (3.1) 

i 

and the LFCC eigenstate is 

ICGP)} = VZe T \<j>:(P)). (3.2) 

We truncate the T operator to just photon emission from an electron 

T= J2 j d y d k±j ' -^V¥t^\yA)aUy^ ± -^b) s {l-y,-k ± -p)bM. (3.3) 

ijlasX 

Surprisingly, this will introduce as much physics as a two-photon truncation of the Fock 
space [19j. The resulting effective Hamiltonian V~ is derived in Appendix IE1 



A. Right-hand eigenvalue problem 



The projection (\2.7\i of the eigenvalue problem onto the valence sector yields, for the 
effective Hamiltonian in flBlTj) . 



2 ± 
m i Z ai 



+ I>i4 = M a 



2 Z ± 



(3.4) 



with a = 0, 1 and M a the ath eigenmass. The self-energy 1^ is defined in fIBlOj) . Clearly, 
the amplitudes z ai are actually independent of spin, and the spin index ± can be dropped. 
We also have left eigenbras of P V V~P V , 



(<^GQI = <o|I>A±(£), 



(3.5) 



for which the amplitudes satisfy 



"lj Zai 



+ Y,(-±) i+i IjiZ«J = M l~Zai- 



(3.6) 



These amplitudes are also independent of spin. 

The left and right valence eigenvectors for different eigenvalues are orthogonal. The 
normalizations are chosen to satisfy 



(C'GP'MGQ) = (-i) a <WW(£' -P) 

Notice that the a = 1 state has negative norm. We then have 

i 

We also have an identity matrix in the valence sector 



'J- 



(3.7) 



(3.8) 



(3.9) 



Projection of the eigenvalue problem onto the one-electron/one-photon sector gives 

1 



£(-!)'** { h ^ x (y> £0 + ^fi\y, k ± ) 



m) + k\ tf x + k\ _ 2 

It VA 



1-2/ 



y 



t^\vM 



To partially diagonalize in flavor, we define 



(3.10) 






c%KvJ±) = £(-i)^^M# a (v,£j-)- 



(3.11) 



With analogous definitions for the vertex functions H and the vertex-loop correction V, and 
with 

(3.12) 



hb> = ( — 1) / jj—lyZbiZb'jIij, 



we have 



Mi 



Ml + fc 2 ^ + fc 2 



1 



(3.13) 



+- 



/«/ 



vs a (»,Sl)-E rr^ A (^ 



Here the eigenmass M& has replaced the bare mass m^ in a natural way, without invocation of 
a sector-dependent parameterization. These equations are to be solved simultaneously with 
the valence-sector equations, f)3.4p and (13. 6p . To facilitate the calculation, the self-energy 
contribution I b y and the vertex correction V° b \ can be expressed directly in terms of the 
wave functions C^ bl x as 



Ihh' 



if £(-i)«*V / ^^(v. ^)^ A (^ £0 (3-14) 

alsX 



and 



y^/fr,^- v ( \Y ,+b ' +l ' e x ' [ d y' d k± Q^-y-y' 

abl ^ ±} " a ,hJ ] J m w-m- 



yf 



TTSS A 



x*f y 



Wa'Xf y 



(3.15) 



y Z* \no<y'\' („i Z.i 



i,i/ r ^_ i"i 



1 + i37*±)C5w (jfr;. ^ + if3*i)^wi' (</, fci), 



where 



H%\y,k±) = J^i-^ZaiZbjK^iy, ki 



K) 



(3.16) 



The original self-energy contribution ly, defined in (1B10[) . can be obtained as 



i i:j = (-iy y^ j (-i) b z bi z b/j i w 

bb' 



(3.17) 



The wave functions C^ x for different J z index a are then seen to be related in the same 
pattern as the vertex functions H° b sX . Since the spin and polarization dependence is not 
affected by the flavor diagonalizations, this pattern can be read from the structure of the 
fundamental vertex functions h a a l x given in (1A22I) . We find 



/~i — hA /->H — A* /~~i A 

U abl ~ °aM ' b aM 



■C: b f X *, for A = ±, 



(3.18) 



^abl 



-C: b J X *, Cfc* = C£* for A = 0,3. 



(3.19) 



B. Left-hand eigenvalue problem 



To obtain the left-hand eigenstates, for use in computation of matrix elements, we also 

need to solve 

— 1.~ .__ M 2 n +P* 



r- ] mp)) 



p^ 



■!<(£)>> 



(3.20) 



with M a fixed and V~ simplified by using £°^ A as a solution to Eq. ( I3.1(jp . for which the 
curly bracket in the third term of (1B17J) is zero. The conjugate is then 



ijs 



V 



5a— — 4 — 

2 _i_ „2 



P 



b] s {p)b is {p) 



(3.21) 



+ J2(-l) l e" f dp_^±^al x (p)a lx (p) 



+ E / dydk L / /-i fi ~3 + h V?( y > fa a \x(Vi k ± ;p)b] s (l - y,-k ± ;p)b itT (p) 



ijlasX ' 



1Gti 3 p 



+ B- iy Ll dvd ^IvA 



dy'dk'j_ 



dp' 



V167T 2 



' p 



'+ 



ijlasX i'l'a'X'" " V iu " "P 

x S((l - y)p + - (1 - y')p' + )8{{\ - y)p x - k ± - (1 - y')p' ± + k' ± ) 
x h%\y, fatf$'*(y', faa] x (y, fc±;p)4^(p')^(p)^A'(^ fag) 

~ E(" 1 ) i E J dydk ± J ^ J dy'dk[hl?(y,k ± )t^'^y\k[) 

ijlasX j'l's'X' 

x oJ A (j/, fc±;p)&j s (l - y, -fap)bji s >(l - y', -fap)a V \>(i/, fap). 

The truncated left-hand eigenvector is 

ICGQ) = IC(Z)) (3-22) 



+ ^|^^^^/^ A (y,^)aJ A (y,^;P)6J s (l- 2/ ,-^;£)|0) ; 



where |0±(P)> = £^ a 4 ± (£)|0> is fixed. 
We diagonalize in flavor, defining 



DZ\y,fa = ^(-i)^^(y,^.), 

i 
•£ = (-1) 6 E(- 1 )^ A f ^^ b t(y,k ± )D^(y,k A 



b'lsX 



and 



E 

a'Vl's'cr'X' 

(-iss'x 1 */ y ti 

XL "bb'l' U > K _L 



ir '+b'+z' e A' / *Mi 0(1 - 2/ - 2/') 



167T3 ^-^(l-y) 



(3.23) 
(3.24) 

(3.25) 



^—faH^{-y— n k ± + -^-,faDz;?(y', fa. 
l-v l-v' l-v' 



This yields 



• 2 _ Ml + k\ _ tf x + fcj 



D&(VJ±) = Hfffakj.) (3.26) 



The right-hand wave functions C^ A are input to this set of equations. The left-hand wave 
functions D°^ x for different J z index are related in the same way as the right-hand wave 
functions: 

D- ah ] X = D+J x *, D- b f = -D+f*, for A = ±, (3.27) 

DJ X = ~D + aU x \ D~ x = D + J x \ for A = 0, 3. (3.28) 

Both sets of equations, left and right, require numerical techniques for their solution. 

IV. ANOMALOUS MAGNETIC MOMENT 

We extract the anomalous magnetic moment of the electron from the spin-flip matrix 
element of the current 

J+ = ^ 7 +^ = 2V+V+ = 2 J2 Mj+- (4- 1 ) 

ij 



The i/ji + are given by (1 A4j) . For our normalization, the general matrix element is [22 



.:>,.*, r, , ,+ ,„,, zz, „.>. ,w »■ , 2, , ^ ± »V 



16^(C(£ + g)|J + (0)|^(P)) = a^W) ± ^^<WW)- ( 4 - 2 ) 



We compute in the Drell-Yan frame 23], where q + = and P /+ = P + , and there are no 
pair contributions, so that 

^ )=2 ?/A/^^- (43) 

We can apply the formalism for the projected expectation value (j2.24p . with the projection 
P s being the projection onto the physical subspace with the two transverse polarizations 
A = 1,2. The current matrix element is given by 

« (P + q)\J + (0)\^(P)) = i \ 1^-^ 44 

- K,WaK - ) ' JdP!(^(P!)\e- T PlPse T \<j>t(P)) 

Note that T is common to both spin projections a — ± and flavor eigenstates a = 0,1; 
therefore, a general expression for off-diagonal matrix elements is not needed here. 

10 



For the chosen truncation, e ±T can be replaced by 1 ± T and only terms up to first order 
in T are to be kept. On substitution of (13.221) for |^), we have 



(r a (P + q)\e- T P}j + (0)P s e T \^(P)) 



167T 3 



$cr±{ z a,0 ~ z al)\ZaO ~ z al) 



(4.5) 



ijj'ls •* K\=± 



J2 eXl at(y^)t^\y,k 



±)Zaj' 



A=0 



for the numerator and 



{^{E)\e- T PtP s e T \<t>t{P)) = 5(P!-P) 



ijls A=0,3 



for the denominator. 

We can then extract the form factors (for a = 0) as 



(4.6) 



1 



F ^-N 



{zqo — ^0l)(^00 — z 0l) 

' dydk ± ] ^ ^ x ^ ^ _ y f ±)f p* {y) %±Uji 



(4.7) 






+/ 



16vr 3 



3 



^%t(yM)^Ky,ki-U 3 



A=0 



and 



Fiiq 2 ) = ±- 



2M 



g 1 ± iq" 



£ ( -i)^-»£(-W# 

jj7 j's 



(4.8) 



u=± 

-E eA ^*(^^)^ A (^^)^ 



A=0 



with 



^ = i - em) ^ e eA / ^s A *^ k)*r (* ^ 



(4.9) 

y7s A=0,3 

However, for opposite spins, / and t are orthogonal, because the azimuthal integration yields 



dydk ± m**(y,k ± )tp*(y,k ± )=0. 



(4.10) 



11 



This eliminates the second term in Fo. In the limit that q 2 — > 0, the first term can be 
rewritten as a derivative, as shown in [22], to obtain a e = ^(0) as 



M, 



".^Ef-D^'+'E 



dydk\ 



ijj'ls A=± ' 

In terms of the wave functions, C^ x and -D^f, we have 

dydkj 
-i) z ai z bj 2_^ J 

ijabls 

and 



d . d 



v l ot(y, *l) orr =F i^) tf,$(y, k ± )z 0j 



dk 1 dk 



(4.11) 



•• = ±# E (-i)*"*- 1 ^ E / w^'^- * 

ijabls A=± 



g , g 



c$£ x (v,ki 



Ms A=0,3 ^ 

In the limit of infinite PV masses, and with Mq = m e the electron mass, 



*i( 



H 



i+E 



71 u=± 

3 

-E eA ^oo A *(^^)4o A (2/^i 



A=0 



and 



with 



^(g 2 



2m e 1 ^-^ ^-^ f dydk\ 



EE 



~q 1 ±iq 2 N ^^ ] 16tt 3 



■^"(y,*±-y?J.)<iw 1 A (ll,£] 



"=i-EE 



a / dydk ± . ± ^ Xi . - .> 



s A=0,3 

In the g 2 — )■ limit, we have -Fi(O) = 1 and 



167T 3 



■&"(.y,ki)t$S(y,k ± ). 



«. = ±?EE/S.«:^f 1 



167T 3 



s A=± ' 

The leading perturbative result is 



d . d 



asXi 



Oy,*x) = Ov.*i.) 



Koo(y:k±) 



2.1.2 „2 , i,2 • 



?/; 



2 mj+fc^ A'pA+fc 



"e 1-y 

When jiq goes to zero, we find J\f = 1 + C(a 2 ) and 

Ot ^, 9 , 

a e = — + 0(a 2 ), 

Z7T 



(4.12) 
(4.13) 



(4.14) 



(4.15) 



(4.16) 



r% r^)4o A Mx). (4.17) 



(4.18) 



(4.19) 



which agrees with the Schwinger result 24j. The result for the normalization factor M 



depends on the fact that h% is the same for A = and 3 when the photon mass \i\ is zero, 



as found in (IA22I) . so that the leading 0(a) contributions to the sum over A cancel. 
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V. SUMMARY 



We have applied the LFCC method [14J to QED in an arbitrary covariant gauge. The 
anomalous moment of the electron is given by Eq. (j4.12p . which must be evaluated with 
use of the left-hand and right-hand functions -D^ A and C^ x . These functions are obtained 
as solutions of the corresponding eigenvalue problems ( I3.26P and (I3.13p . The leading per- 
turbative contribution to the anomalous moment is verified to be the standard Schwinger 
term; a nonperturbative solution requires numerical methods. To carry out this construc- 
tion, we have extended our analysis of arbitrary gauges |9J to include a specific choice of 
projection onto physical states, as discussed at the end of Appendix |Aj and have extended 
the calculation of LFCC matrix elements jl_4j to include such projections, as expressed in 
Eq. (H25). 

This provides an extensive test of the LFCC method in a gauge theory regulated by PV 
fields. However, the method is much more general than this; it should be applicable to any 
regulated light-front Hamiltonian. The avoidance of a Fock-space truncation then provides 
several benefits, with the absence of Fock-sector dependence and spectator dependence and 
of the uncanceled divergences that can result from them. The method allows for systematic 
improvement of a calculation, through the addition of terms to the exponentiated operator 
T. Thus a natural next step is to include one or more positron terms, to study the dressed- 
photon state [2l|, pair contributions to the dressed-electron state, and positronium. 
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Appendix A: Light-front QED in an arbitrary gauge 

Here we summarize a formulation of light-front QED regulated by the inclusion of Pauli- 
Villars (PV) electrons and photons [26] and quantized in an arbitrary covariant gauge. 
Details can be found in [9|. We write the Lagrangian as 



1 """ " I,, 2 4^ 4. _ I, 



- 7 rV + nA^A,, ~ n( (3*A„) 2 



£ = £(-! 

i=0 
2 



(Al) 



j=0 



with 

2 2 

V' = X)V / AV'i, ^ = J2V& A ^ Fi^ = a^Ah, - dvAip. ( A2 ) 

The index % is zero for physical fields and 1 or 2 for PV fields. The coupling coefficients $ 
and £j are constrained to satisfy % — 1, £0 = 1, and 

2 2 

Ec- 1 )^ - E(-w = °> ( A3 ) 

i=0 i=0 
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and to reproduce the correct chiral symmetry in the limit of a massless electron [8] , and to 
guarantee a massless eigenstate for the photon [21] . 

The fermion fields ipi are decomposed into dynamical and nondynamical parts ipi± = A±ipi 



o„,± 



by the complementary projections A± = 7^7 

1 



/2 [3|, |25]. The dynamical part is written 



A 



VlQ 



7T" 



J2 / rf ^ 6 fa (&)e-* ,2 + <-.(fc)e** 



(A4) 



with 



X- 



1 

71 



/1\ 



1 
W 



X- 



1 



/ \ 

1 



\- l J 



and 



R.(*0,d!v(*'} = (-1)%'<W%- £')■ 
The nondynamical part satisfies the constraint 



(A5) 

(A6) 

(A7) 

(A8) 



(27V 



-1)^^+ -ie.A L < v /&$^V'i- 



•1)^7°^+. 



However, the constraint for the sum that enters the interaction becomes simply 



(A9) 



where the photon field does not appear. The constraint can then be trivially solved and the 
nondynamical field removed from the Lagrangian. 
A vector field of mass \i\ is written as 



Ai„(x) 



dk 



VL&nHr 



J2^\k)\a lx (k)e- ik - x + al(k)e ik - 3 



(A10) 



. A=l 



+ef{k) \a l0 (k)e- fkx + a| (fc)e 4 ^ 

with k a four-vector associated with a different mass jxi = fii/y/C, such that 

~k = k, k~ = {k 2 ± + ti)/k + . 



(All) 



This allows the field to satisfy the Euler-Lagrange equation, while the fourth polarization 
A = does not satisfy the gauge condition d-A[ = J9|]. The polarization vectors are defined 
by 



e (1 ' 2) (fe) = (0,2ei, 2 .^j./A;+ei l2 ), 

e^(k) = ((kl- tf)/k + ,k + ,k ± )/» h 
e {0 \k) = klin = ((k 2 ± + tf)/k + , k+, k ± )/H, 



(A12) 
(A13) 
(A14) 
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and the commutation relations are 

[aix(k), 4 A ,(fc')] = (-l) l Sii'e x S X y6(k - A;'), (A15) 

with e=(— 1,1,1, 1) the metric signature for the physical photon. 

For the purpose of making J 2 -conservation explicit, it is convenient to introduce circular 
polarizations 

e^ = T^{e {l) ±ie^) (A16) 

and the associated creation and annihilation operators 

a \± = T-^= OJi ± ia\ 2 ), a l± = =F-/=(an T ia a ). (A17) 

These operators satisfy the same commutation relations (IA15I) . with A now taking the values 
± instead of 1 and 2, and with e = 1. Sums over A will, in general, include 0, ±, and 3. 

The light-front Hamiltonian is then given by 

V-=V 0a + V ob + Vr nV (A18) 

with 

-P 0a = E(" 1 ) /eA f dp^^al(p)a lx {p), (A19) 

ix J P 

-Pot, = E(" 1 ) i f dt^^bl(p)b is (p), (A20) 

is •* P 

and, if antifermion terms are neglected, 

T M = J2 J d y d k J -M= (A21) 



1 f)7T T) 
ijlasX u V " 

x \Kji X (y> k±)4 x (y, k x ;p)b] s (l - y, -k±;p)b i(r (p) 
+h™ x *(y, k ± )bi(p)b js (l - y, -kr,p)aa{y, k±;p) 



As discussed at the beginning of Sec. El the multiple arguments of the creation and anni- 
hilation operators are defined by a,i\(y,k±;p) = ai\(yp + ,yp± + k±). The vertex functions 
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arc 



-> k i p^^ 

hfi±{v,k ± ) = Te^mMi (l ^ )3/2 , (A22) 



y)z/ 3/2 

■ k ± e ±i<t} 

y 

rrii(l — y) — m 



hf^ivA) = ±e^/2PMi—T 






±±3/. z \ „ at^t- ™^y 2 - /^ 2 (! - y) + *i 



/J-z(l — 2/)2/ 3/2 

. ±±o r?/ r n _ p /Tr^ m,m^ 2 + g(l - y) + fcj 

m l -y)vy 

The gauge condition <9 • A; = for the Ith flavor is implemented as a projection onto a 
physical subspace. Only the polarization with opposite metric, A = 0, contributes to d ■ Af, 
we have 

. 4 = -,'i| | -^M= [a l0 (k)e- fk -* - aj^e**] . (A23) 

We require that physical states iV'phys) satisfy 

(^ phys \d ■ A t \tp phys ) = 0, (A24) 

which is guaranteed if |?/>phys) is annihilated by the positive-frequency part of d ■ A t or, 
equivalently, 

Hia l0 (k)\ip phys ) = 0. (A25) 

However, in the massless limit, the polarization vectors e^ and e^ become identical. This 
suggests that implementation of the gauge projection for a massless photon should involve 
the removal of not just the e^ polarization but also the e® polarization, as one would 
expect on physical grounds. 

To make this more precise, define two new polarizations 

g(0,3)(£) = _^L (e (3) + e (0 }) ± ^ (e (0) _ e ( 3)) (A26) 

2k^ 2fii 

before taking the massless limit. With use of the definitions in Eqs. ( 1A13I) and (IA14I) . they 
can be written explicitly as 

- ( „,s, (i) = !_ ^l + ft '(l-0/2C±(^(l + C )/2C tt j (A2?) 

These have the useful properties of not being the same in the massless limit and of reducing, 
in Feynman gauge (C = 1), to the standard choice of [8| 

_ (3) = k-(n-k)n ^ g(0) = 

n ■ k 
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where n is the timelike four-vector that reduces to n — (1, 0, 0, 0) in the frame where k± = 0. 
By requiring that 

J2z (X) (k)aix(k)= J> (a) (£Ha(£), (A29) 

A=0,3 A=0,3 



we find that 

and that the a\ X satisfy the same commutation relations as do the a\\. 
The gauge projection ( 1A25I) now becomes, in the massless limit, 



(A30) 



{aio(k) - a l3 (k))\i; phys ) = 0. 



(A31) 



This removes from \ip p hys) one null combination of aio and a/3. The other null combination, 
a^o + S/3 remains but makes no contribution to physical quantities simply because it is 
null; instead, it represents the remaining gauge freedom not fixed by the Lorentz gauge 
condition |8j. For practical calculations of observables, it too can be removed, which makes 
the gauge projection equivalent to the removal of both the e^ and e^ 3 -* polarizations. 

For the massive PV photons, we apply the same projection. This could not be done 
for a physical massive vector particle, for which three polarizations must be retained. For 
the PV photons, which are regulators and not physical, the projection becomes part of the 
regularization prescription. In any case, in the infinite-mass limit, all such contributions 
from the PV photons disappear; removing them prior to taking this limit will not change 
the result. 



Appendix B: Derivation of the Effective Hamiltonian 

The effective Hamiltonian V~ is constructed with use of the Baker-Hausdorff expansion 
( 12. 9p as applied to V~ in (IA18B and T in (13.31) . The first commutator with the first term in 
V~ is 



[-Poa,T]= J2 [ d v&±[-A=. 



Ma + (yp± + ^J 



x£^ A (2/, k±)a] x (y, k ± ;p)b} 3 (l - y, -k±;p)b il7 (p). 
For the second term, we have 



i^HE/^/^ 



m) + ((1 - y)p± - k ± f 

^y 



{m 2 i + p\) 



(Bl) 
(B2) 

(B3) 



xt™ t x (y, k±)a] x (y, k±;p)b] s {l - y, -k±;p)b i(T {p) 



Due to cancellations between terms involving the transverse momentum p±, the combination 
of the two simplifies to 



\^a + V ob ,T}= J2 [dydk ± f 



dp 



ijlasX ' 



a/167T 3 ^ 



m} + k 2 ± | tf x + k 2 ± 


2 




- m\ 


i-v y 





(B4) 



xt-j a i X (y, k±)aj x (y, k ± ;p)b\ s (l - y, -k±;p)b ilT (p). 



17 



A graphical representation is given in Fig. [TJ Since this commutator already involves a 
net increase by one particle, no additional commutators with T need to be considered for 
these terms, in the chosen truncation. Similarly, the photon emission term in pr t already 
increases the particle number and need not be considered at all in any commutator. For the 
photon absorption term, however, two commutators must be considered. 



/ \AXW 





(a) 



-«— x — «- 





-«— X 




< X < 



(b) 

FIG. 1. Graphical representations of the terms in (a) [Po~ a ,T] and (b) [P^T]. Each diagram 
represents an operator that annihilates particles on the right and creates particles on the left. The 
crosses represent kinetic-energy contributions. 

The first commutator with V^ t generates terms that do not change particle number. 
They can be represented graphically as in Fig. |2j The commutator can be written as 



-* — *=) — * — ' 





* h <- 



71 + 71 + 75 + 71 



71 




71 



•*-? * *- *- 




71 



<-<= — « — ■=>-*■ 



r 4 




< — H — zr* b — *- 




FIG. 2. Graphical representation of the terms in [7-' i t , T^| 



[V^T} = T 1 + T 2 + % + %, 



(B5) 
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with 






dp' 



v^ 



// 



'+ 



7T' ; 



ijlasX i'j'a's' V " 

x htf*{y, k±njn\y', k' ± )6(yp + - y'p' + )5(yp ± + k ± - y'p' ± - k' ± 

x b), s ,(l - y', -k^j/)b\ a (p)b js (l - y, -£±;p)W(p'), 
dp 



(B6) 



t 2 = £(-W' / "Mm^ApI 

ijira' 

Ts = - E^" 1 ) 1 E j dydk^ J dy'dk' L j ^h^\yM)t^\y\k' L ) 

ijlcrsX j'l's'X' 

x 4y(V> JaJlO&JvC 1 ~ 2/'> -ft^M 1 - V: -k±;p)aix(y, k±;p) 

* - EW / ** / Tifc E / <v* / t! 



/p, 



tf«A •' - \/ 167r V , — v J J V107T- 

x <5((1 - y)j9 + - (1 - y')p' + )8((l - y)^ ± - fcj. - (1 - y')i? ± + k' ± ) 
x h™ x *(y, k ± )t^'(y', k[_)4 x ,(y', fc|;p>J ff (p>^(p> A (y, kr,p), 



(B7) 

(B8) 

(B9) 



and 



^ = (-i)*E(- 1 ) i ' +,eA / 

J'/sA "^ 



^^f(y,^)4f(y,^). 



(BIO) 



The first term, 7i, does not contribute to the chosen truncation to one fermion. The spin 
dependence of the self-energy term T2 can be simplified upon the observation that the 
dependence of t^/f A (y, k±) on the azimuthal angle of k±_ must follow the pattern of dependence 
in h^/^y, k±), in order that T conserve J z . The azimuthal integral in IJ° then implies that 



jaa e j 

-Lji O a(T >lji, 



(Bll) 



with I™ real and independent of the spin projection a. This then simplifies T2 to be 



ijrs J ^ 



(p)hs(p)- 



(B12) 



The second commutator [[P^ t ,T],T] generates several more terms, all of which increase 
particle number by one; however, most of these do not contribute. Graphical representations 
are given in Fig. |3j In particular, [7i,T] contributes nothing for our truncation. The next 
possibility [%, T] does contribute fully and, of course, has a structure very similar to that 
oi[V Qb ,T]: 



ijlasX ' 






16ttV 

x a} x (y } k ± ;p)b] s (l - y, -k ± ;p)b irJ (p). 



(B13) 
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* — C? «-? < r-* 





* — 7 — * b— < s— * 




(c) 

FIG. 3. Graphical representation of the contributions from [['P^ t ,T , ],T], with (a), (b), and (c) 
corresponding to [72, T], [T3, T], and [74, T], respectively. [7i,T] does not contribute at all for the 
chosen truncation. 



The third term contains several pieces but contributes only 

dp 

I ayan ± 1 — 






167T 3 p H 

x aJ A (y,fcj.;p)6j s (l -y, -k±;p)b ia (p), 



(B14) 



which combines with the second term of [72, T]. The fourth also contains several pieces; all 
that remains in the given truncation is 

[T4,Jl-> £ [ dydk ± [ ^^=V^Hy,k ± )al(y,k ±] p)b] s a-y^-k±;p)bUp), (B15) 

.rrrw y \/i67r%+ 



ijlusX ' 



where 



I/<" A f7,^- V f iV'+J'+V f dy^hi d{l-y-y') 



(B16) 



1 — 1/ J 1 — 1/' 1 — 1/' 



i-y' 



This is a vertex correction, as illustrated in Fig. [3]^c). 

The effective Hamiltonian for the given truncation is then the sum of all these terms, 
according to the Baker-Hausdorff expansion. A graphical representation is shown in Fig. HI 
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The final expression is 



b] s (p)b is (p) 



(B17) 



+ J2(-l) l e x f dp 
ix •* 

+ 2_\ / dydk± / 



Phx+pIj 



p 



<H\(p)a>i\(p) 



dp 



1, 



a/167T 3 P 

+ 



h^(y,k ± ) + ^\y,k ± ) 



m^ + kl , ^ + A;i 



i-y 



+ 



- to, 



*«!*(!/, £l) 



+ \^ / dydk± / ■ 



4 E tz^*i(v> £o - D- 1 )*'"^ Sow 

x a\ x (y,k ± ;p)b} s {^-y,-k±;p)b ia (p) 

dp , -. . 

=^Kji *(j/> fcjJ&lrUyM 1 - 2/, -k±;p)aix(y, k±;p) 



ijlsaX 



dp 



lQ n 3p+ 



dy'dk\ 



dp' 



V 167T 3 



'P' 



16"7T 3 p 

+ E(- 1 )' E j dydk^j 

ijlasX i'l'a'X' 

x tf((l - y)p + - (1 - yV + )<S((l - v)p± - fc± - (1 - 2/')p1 + ft) 
x ^f*(y, EOtf/d/, fci)4 v (2/', £i;p>l ff (#^(P>A(y, £i_;p) 

^ E^" 1 ) 4 E /^ /i^s / dy'dk' ± h^ x *(y,k ± K S '7(y',k' ± ) 

ijlasX j'l's'X' J J J 

x 4a'(s/ / >&;p) 6 5v( 1 -lA-fe^M 1 -y,-k±;p)aix(y,k±;p) 
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FIG. 4. Graphical representation of the effective Hamiltonian V , limited to those terms that 
contribute in the given truncation of the operator T. 
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